function[separacion]=cohesion(xbin)
[m,n]=size(xbin);
xx=zeros(m,n);
xx(find(xbin==0))=1;
[L,num]=bwlabel(xx,8);
separacion=[];
distancia=[];
anterior=0;
a=0;
if num>1
    for s=1:num
       [f,c]=find(L==s);
       if (length(f)>10)&&(a==0)
           anterior=max(c);
           a=1;
       elseif (length(f)>10)&&(a~=0)&&(min(c)>anterior)
           distancia=[distancia (min(c)-anterior)];
           anterior=max(c);
       end
    end
end
if isempty(distancia)==0
    separacion=mean(distancia);
end